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ABSTRACT 

We combine new Spitzer Space Telescope observations in the mid- and far-infrared with SCUBA 
850/im observations to improve the measurement of dust temperatures, masses and luminosities for 11 
galaxies of the SCUBA Local Universe Galaxy Survey (SLUGS). By fitting dust models we measure 
typical dust masses of 10 ''■^M0 and dust luminosities of '--^IO^^Lq, for galaxies with modest star 
formation rates. The data presented in this paper combined with previous observations show that 
cold dust is present in all types of spiral galaxies and is a major contributor to their total luminosity. 
Because of the lower dust temperature of the SCUBA sources measured in this paper, they have flatter 
Far-IR i^Fj^(160/im)/^F^(850/j,m) slopes than the larger Spitzer Nearby Galaxies Survey (SINGS), the 
sample that provides the best measurements of the dust properties of galaxies in the nearby universe. 
The new data presented here added to SINGS extend the parameter space that is well covered by local 
galaxies, providing a comprehensive set of templates that can be used to interpret the observations of 
nearby and distant galaxies. 

Subject headings: Galaxies: fundamental parameters - galaxies: ISM ~ galaxies: photometry 



1. INTRODUCTION 

The mid and far infrared 10-1000 fim) is a wave- 
length domain where some galaxy populations show sig- 
nificant changes over time (e.g., Puget et al. 1996; 
Sanders & Mirabel 1996; Chary & Elbaz 2001; Elbaz 
et al. 2002; Le Floc'h et al. 2005; Lagache, Dole & 
Puget 2005) and is therefore critical in our understand- 
ing of how galaxies form and evolve. The strong emission 
of galaxies in these wavelengths has been measured for 
objects with very different optical properties (spirals, el- 
lipticals, Seyferts, interacting galaxies). The optical to 
IR ratio varies strongly from galaxy to galaxy and can 
be as low as few percent for the extreme IR galaxies 
{Lifj > IO^^Lq) where the IR emission dominates their 
total luminosity (Rieke & Lebofsky 1978). 

This IR emission is mostly due to starlight absorbed by 
dust grains and re-emitted at longer wavelengths, though 
some objects have a significant non-thermal contribu- 
tion coming from ionized gas in the interstellar medium 
(ISM) (Laurent et al. 2000a) and from active galactic 
nuclei and radio sources. The mid-IR light (A 5-30 
/im) is emitted by dust grains of different sizes (Desert, 
Boulanger & Puget 1990), that include very small grains 
(VSG), which are responsible for the thermal emission 
at 10-16 /im and are located in HII regions (Laurent et 
al. 2000b), and aromatic molecules (e.g., PAH, Leger 
& Puget 1984) that absorb ultra-violet photons gener- 
ated in photo-dissociation regions (PDR) and re-emit in 
several bands between 3 and 17 /zm (e.g., Laurent et 
al. 2000b; Lagache, Puget & Dole 2005). At longer 
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wavelengths the emission is due to larger grains com- 
posed mainly of silicates (Desert et al. 1990). The latter 
absorb light from the ISM that is generated by young 
OB associations and from the old disk population (Lons- 
dale & Helou 1987). The ubiquity of dust in galaxies 
was demonstrated by the Infrared Astronomical Satel- 
lite {IRAS) survey, which covered over 90% of the sky in 
wavelengths ranging from 12 /im to 100 /im, and showed 
that most of the detected infrared luminosity could be 
explained by dust at temperatures of 40K and higher. 
However, evidence for dust at more than one temper- 
ature was first noted by Helou (1986) from the IRAS 
12/25 /im to 60/100 /im color-color diagram which sug- 
gested that in addition to dust excited by the actively 
star forming regions, another component due to cooler 
"cirrus" was present and motivated the models proposed 
by Lonsdale & Helou (1987) mentioned above, and val- 
idated by Sauvage & Thuan (1992) in study comparing 
the Ha and far-IR luminosities. From the analysis of the 
gas to dust ratio in spiral galaxies, Devereux & Young 
(1990) could explain the observed discrepancy between 
the gas to dust ratios measured in external galaxies with 
the measurements obtained for the Galaxy by assuming 
the presence of a diffuse and cold (T < 30K) dust compo- 
nent dominant in the outer regions of galaxies, while the 
warmer dust would be predominant in the inner regions 
of the galaxy. 

The detection of this cold dust required reaching longer 
wavelengths than accessible with IRAS. This became 
possible with the flight of the Infrared Space Observa- 
tory (ISO, Kessler et al. 1996) which allowed probing 
out to 200 /im and with the commissioning of the Sub- 
millimeter High Angular Resolution Camera (SHARC) 
at the Caltech Submilliter Observatory (e.g., Serabyn 
et al. 1999) and on the James Clerk Maxwell Tele- 
scope of the Submillimetre Common-User Bolometer Ar- 
ray (SCUBA) (Holland et al. 1999), observing at 850 
/im. 

The ISO observations of the nearby galaxies M 31 
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TABLE 1 
Galaxy sample 







T 




T no- f J \ 


T no- ^^ \ 


SFR 










h . 


h . 


A.L . \ 1 


IC 797 


0.00320 


6 


-16.81 


9.27 


8.81 


0.2 


IC 800 


0.00779 


5 


-18.83 


9.07 


9.38 


0.4 


NGC 803 


0.00703 


5 


-19.28 


9.37 


9.63 


4.5 


NGC 3270 


0.02112 


3 


-21.22 


10.23 


10.56 


2.0 


NGC 4712 


0.01473 


4 


-20.91 


9.87 


10.15 




NGC 6137 


0.03153 


-5 


-21.98 


<9.94 


<9.54 


2.9 


NGC 7047 


0.01957 


3 


-20.50 


9.98 


10.35 


3.5 


NGC 7653 


0.01433 


3 


-20.48 


10.17 


10.43 


0.7 


PGC 035952 


0.01331 


4 


-19.73 


9.61 


9.80 


1.0 


UGC 10205 


0.02211 


1 


-20.80 


10.10 


10.35 


2.9 


UGC 12519 


0.01471 


5 


-19.76 


9.95 


10.30 


2.7 



Note. — Redshifts and morphologies come from the NASA Extragalactic Database (NED). The absolute magnitudes were derived from B magnitudes 
in NED. The Lj-r measurements come from Vlahakis et al. (2005), while the Ltjr are derived from the the fits described in Section 3. The star 
formation rates are estimated following Bell (2003) and have estimated uncertainties ~ 50%. 



(Haas ct al. 1998), M 33 (Hippclcin ct al. 2003) and 
M 101 (Tuffs k Gabriel 2003; Sauvage, Tuffs & Popescu 
2005) and by Serabyn et al. (1999) for NGC 891 using 
SHARC, in wavelengths above 100 nm provided the first 
direct detections of cold dust. In a study of NGC 891, 
Popescu & Tuffs (2003) showed that the cold dust ex- 
tends beyond the stellar disk, a similar behavior observed 
in HI disks, thus confirming the scenario proposed by 
Devereux & Young (1990). The data acquired for M 31 
(Haas et al. 1998), M 101 (Tuffs & Gabriel 2003), M 33 
(Haas ct al. 1998; Hinz ct al. 2004) and NGC 55 (En- 
gelbracht ct al. 2004) confirm that the cold and warm 
dust arc distributed differently throughout the galaxy. 
In general, the warm dust component is easily identified 
in spiral arms and nuclei thanks to the larger density 
of star-forming regions, while the cold component is the 
source of the diffuse emission in the disk, though some 
localized components exist in the spiral arms (Tuffs & 
Popescu 2005). 

The ground based observations using SCUBA to target 
nearby galaxies with IRAS detections also suggested the 
presence of significant amounts of cold dust. Dunne et 
al. (2000) obtained SCUBA observations for 104 galaxies 
and were able to obtain measurements of the local sub- 
millimeter and dust mass functions. However, because 
the original sample of Dunne ct al. (2000) was selected 
on the basis of IRAS fiux, they would preferentially de- 
tect galaxies with a significant amount of dust at T>40K. 
Based on the high gas-to-dust ratio (~ 600) compared to 
the Galactic value, Dunne et al. (2000) suggested that 
a cold dust component would be present in their sam- 
ple, but which would not be directly detected by IRAS. 
To overcome the bias of using only /i?AS'-sclcctcd galax- 
ies, Vlahakis, Dunne & Bales (2005) used the CfAl red- 
shift survey (Huchra et al. 1983) to define a sample to 
observe with SCUBA based exclusively on their optical 
properties. These were observed at 850 fxm as part of 
the Scuba Local Universe Galaxy Survey (SLUGS). The 
SLUGS sample confirmed previous findings (e.g.. Tuffs 
& Popescu 2003) using ISO that most of the dust mass 
is contained in cold particles. Whenever SCUBA data at 
450 /im were available, supplementing the 850 /im and 
IRAS measurements. Vlahakis et al. (2005) found that 
two dust components (one warm, another cold) gener- 



ally produced better fits to the observed distribution of 
Mid- to Far-IR fiuxes. Their fits suggested that the total 
cold dust mass in these galaxies could be up to 100 times 
larger than the dust responsible for the emission at 5 /um 
through 70 /zm. 

Our understanding of dust properties of galaxies in 
the nearby universe has increased dramatically with the 
Spitzer Nearby Galaxies Survey (SINGS) key project 
(Kennicutt et al. 2003). SINGS has obtained exten- 
sive observations ranging from the far-UV observed with 
the Galaxy Evolution Explorer (GALEX, Martin ct al. 
2005) to the far-infrared for a representative sample of 
nearby galaxies comprising a wide range of luminosities, 
sizes and morphologies (e.g.. Dale et al. 2005; Dale et 
al. 2007; Draine et al. 2007; Smith et al. 2007). These 
objccits were pre-selected considering a parameter space 
defined by the optical morphology, the luminosity and 
the optical/FIR luminosity ratio (Kennicutt et al. 2003). 
Because prior to the launch of Spitzer most of the avail- 
able mid- and far-IR measurements used to define the 
sample were obtained by IRAS, mainly sensitive to the 
emission by hot and warm dust, it is possible that galax- 
ies with large amounts of cold dust are under-represented 
in SINGS. To test if this bias is present, we will use 
galaxies from SLUGS to address the far-IR properties of 
galaxies with a large cold dust component, and to pro- 
vide a more complete view of the properties of galaxies 
over the whole extent of the IR wavelengths. In Section 
2 we present the sample and describe the observations. 
In Section 3 we use modified black-body curves and dust 
models to estimate the temperatures of the dust compo- 
nents and characterize the Spectral Energy Distributions 
of these galaxies. Section 4 presents our conclusions. An 
Ho= 70 kms~^Mpc~^ , Qm = 0.3 and Oa=0.7 cosmology 
is assumed throughout. 

2. DATA 

2.1. Sample definition 

The Spitzer observations targeted a subsample of 11 
galaxies from Vlahakis et al. (2005) that have SCUBA 
observations in 850 /zm (and in three cases additional 
data at 450 ^m). The sample, presented in Table 1, con- 
tains galaxies with morphologies ranging from elliptical 
to late spiral, but with no cases of irregular or interacting 
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TABLE 2 

Detector Parameters 



Detector 




Pixel Size 


FWHM 


Conversion Factor 


Calibration Uncertainty 




fini 


arc sec 


arc sec 


Jy/(MJy/Sr) 




IRAC 3.6 


3.550 


0.86 


1.66 


1.7492x10-5 


1.8% 


IRAC 4.5 


4.493 


0.86 


1.72 


1.7492x10-5 


1.9% 


IRAC 5.8 


5.731 


0.86 


1.88 


1.7492x10-5 


2.0% 


IRAC 8.0 


7.872 


0.86 


1.98 


1.7492x10-5 


2.1% 


MIPS 24 


23.68 


2.49 


6.0 


6.514x10-6 


4% 


MIPS 70 


71.42 


9.85 


18 


1.60 


5% 


MIPS 160 


155.90 


16 


40 


0.25 


15% 



galaxies. All galaxies in this sample are relatively quies- 
cent in terms of their star formation activity with ~ 2 Mq 
yr~^, estimated using the Bell (2003) relation. Because 
these galaxies are at moderately low redshift, with typ- 
ical sizes of ^ 2', the main subcomponents (bulge, disk, 
possibly spiral arms) are easily resolved in data ranging 
from the optical down to 24 ^m. 

2.2. Observations 

The observations were taken as part of the Spitzer 
Guaranteed Time Observation Program 30348 (P.I.: G. 
H. Rieke) and used both the Infrared Array Camera 
(IRAC, Fazio et al. 2004), which is sensitive to the 
Rayleigh- Jeans stellar and warm dust emission, and the 
Mid Infrared Photometer for Spitzer (MIPS, Rieke et 
al. 2004), which is sensitive to warm dust in its bluer 
wavelength windows (24 /im and 70 /im) and cold dust 
at 160 /im. 

The IRAC data were taken using all four wavelengths 
(3.6, 4.5, 5.8 and 8.0 /im), where a 5' x 5' region of the 
sky was observed, using a 5 position dithered pattern, 
each exposed for 30 seconds, giving a total exposure of 
150 seconds per band. 

The MIPS observations were taken in photometry 
mode, where the spacecraft points at the object and then 
performs a series of mirror motions and spacecraft ma- 
neuvers to dither the source over the array (Rieke et al. 
2004). All MIPS 24 /im observations were taken in a 
single cycle using a sequence of 12 dithers, followed by 
another series of 12 images which are offset to calculate 
the sky background. The total on-source integration is 
120 seconds. Most 70 fj,m observations were also taken in 
a single cycle, though for two galaxies, because of their 
large apparent extent, "fixed cluster" observations using 
two cycles were necessary; the total integration time per 
pixel amounts to ~ 125 seconds. In the case of 160 /^m, 
the observations were taken over 4 cycles, which adds to 
a total integration time of ~ 42 seconds per pixel. 

Table 2 summarizes the instrumental parameters for 
each detector, consisting of the nominal wavelength (e.g.. 
Reach et al. 2005; Rieke et al. 2008), the pixel size in the 
final mosaic, the FWHM of the point-spread function, 
the factor used to convert fluxes into Janskys and the 
calibration uncertainty. 

2.2.1. Data Reduction 

The pipeline-processed IRAC Basic Calibrated Data 
(BCD) were retrieved from the Spitzer Science Center 
(SSC) and were post-processed using the IRACproc suite 
of scripts (Schuster et al. 2006). These combine the SSC- 



processed images, reject outliers, correct for bad-pixels 
and cosmic rays, and produce the final flux-calibrated 
mosaic (Schuster et al. 2006), using the procedure and 
calibrations of Reach et al. (2005). The flnal photometry 
includes the calibration corrections determined by Rieke 
et al. (2008). The photometric uncertainties in the IRAC 
calibration are 1.8%, 1.9%, 2.0% and 2.1% for 3.6 /im, 4.5 
/im, 5.8 fim and 8.0 /im respectively (Reach et al. 2005). 
The IRAC measurements already contain the corrections 
for extended sources^. No color corrections were applied 
to the tabulated data. 

The MIPS data were reduced using the Data Analysis 
Tool version 3.06 (DAT) which is extensively described 
by Gordon et al. (2005). The DAT reads the raw data 
downloaded from the SSC archive, converts the integra- 
tion ramps into slopes (in data numbers per second) re- 
moving at the same time instrumental effects, electronic 
non-linearities and stimulator flash latents in the case of 
70 and 160 /im observations. The DAT also performs 
dark current, flat field and illumination corrections, re- 
moves distortions introduced by the off-axis location of 
the detectors, and creates a flux-calibrated mosaic (Gor- 
don et al. 2005). The calibration of the MIPS data, 
described by Engelbracht et al. (2007), Gordon et al. 
(2007) and Stansberry et al. (2007), used repeated obser- 
vations of stars which are referred to a model spectral en- 
ergy distribution of an ideal A star ("Vega"), calibrated 
using the values in Rieke et al. (2008). The calibration 
uncertainties are ^ 4% for 24 /im (Engelbracht et al. 
2007), 5 % for 70 pm. (Gordon et al. 2007) and - 15 % 
for 160 /im (Stansberry et al. 2007). 

2.2.2. Ptiotometry 

Galaxy fluxes were measured using the IRAF^ ellipse 
task (Busko 1996), based on the method of Jedrzejewski 
(1987) to fit elliptical and circular isophotes. This task 
takes first guess values of object center, semi-major axis, 
position angle and ellipticity to calculate the ellipse at 
the starting radius. Excepting for the semi-major axis, 
the fitted parameters are used as input values for the 
next radius (outer or inner), which is calculated using 
a radial step set by the user. The ellipse task allows 
masking bad pixels as well as other objects in the image, 
such as foreground stars and background galaxies. Since 
ellipse does not calculate sky background values, these 

^ |http : // SSC ■ spitzer ■ caltech.edu/lrac/callb/extcal /index . html] 

''' 1HA1<' is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA), Inc., under cooperative agree- 
ment with the National Science Foundation. 




Fig. 1. — Postage stamp images of galaxies in this sample. The galaxy identification is noted at top of window, while the band is noted 
in the lower right corner. A scale bar corresponding to 2 kpc at the galaxy distance is shown in the upper left corner of each figure. 



















PGC036S5: 




PBCCI35952 








PGC 03696: 












g 












g 














i 

e 

¥ 


«> 

. .1 .... 1 .... 1. . 


! 

¥ 


5.8 ;im . 

..).... 1 .... M . r 


¥ 


• 

B.D . 




2+ ^im , 
. .1 .... 1 .... 1 . . 


1 

e 

? 


70 ^ . 
. . 1 . . . . 1 . . . . 1 . . 


¥ 


^ ^ \ 

ISO ^tfn . 

. . 1 . . . . 1 . . . . 1 . . 






» -50 




50 -so 




50 -M 




» □ -50 




50 -50 




50 -50 




50 □ -50 






IJX1 0205 




UGC10205 




UCC10205 




LtGC102iJ5 




UCC1 0205 




UGC 10205 




UCC1020S 












g 
























• # 




m 


I- 












V 


m 










i.e /xn . 


¥ 


4.5 . 


? 


S.S _ 


¥ 




¥ 


24 ^rr . 




70 ^im . 


¥ 


I6D fun . 






so a -so 
4ft Me 




50 -» 




50 -50 




40 □ -H 




50 -50 
«e Ml 




50 -ao 

Vt AH 




sa □ -50 






U0C12519 




UGC12519 




UGC12519 




IJSC1 2519 




UGC12519 








LGC12519 
































1° 




/ 




/ 




f 




/ 






















4.5 . 


s 

1 




1 


6,0 ^iin _ 


¥ 




£ 


70 . 


il 
1 


190 fan . 






w -so 








H -BO 

are *K 




U -M 
are w 












so -eg 



5 



Fig. 1. — (cont.) 



TABLE 3 

Galaxy sample: Spitzer Measurements 



Name 




3.6/jm 




4.5/im 




5.8/im 


8.0Aim 


24/jm 






70 m 


m 






160/^m 






Jy 






Jy 






Jy 




Jy 




Jy 






Jy 






Jy 


IC 797 





021± 


003 





014± 


002 





027± 


004 


0.065± 


010 


0.072± 


007 





939± 





188 


2 


708± 0.677 


IC 800 





017± 


003 





012± 


002 





020± 


003 


0.042± 


006 


0.046± 


005 





519± 





104 


1 


463± 0.366 


NGC 803 





045± 


007 





029± 


004 





042± 


006 


0.092± 


014 


0.071± 


007 





904± 





181 


4 


241± 1.060 


NGC 3270 





050± 


007 





031± 


005 





043± 


006 


0.103± 


015 


0.083± 


008 





725± 





145 


3 


727± 0.932 


NGC 4712 





028± 


004 





019± 


003 





027± 


004 


0.071± 


Oil 


0.069± 


007 





711± 





142 


3 


236± 0.809 


NGC 6137 





058± 


009 





036± 


005 





022± 


003 


0.018± 


003 


0.007± 


001 





017± 





003 





030± 0.007 


NGC 7047 





028± 


004 





019± 


003 





031± 


005 


0.067± 


010 


0.060± 


006 





722± 





144 


2 


201± 0.550 


NGC 7653 





046± 


007 





030± 


005 





060± 


009 


0.145± 


022 


0.131± 


013 


1 


776± 





355 


4 


346± 1.087 


PGC 035952 





011± 


002 





006± 


001 





Ollib 


002 


0.028± 


004 


0.045± 


005 





611± 





122 


1 


211± 0.303 


UGC 10205 





042± 


006 





028± 


004 





029± 


004 


0.056± 


008 


0.041± 


004 





727± 





145 


1 


762± 0.441 


UGC 12519 





020± 


003 





013± 


002 





033± 


005 


0.078± 


012 


0.108± 


Oil 


1 


081± 





216 


3 


690± 0.922 



TABLE 4 

Galaxy sample: Measures from the Literature 



Name 




J 

mag 






H 

mag 






K 

mag 






12/jm 

Jy 


60/im 

Jy 






100/jm 

Jy 


450/^m 

Jy 




850/jm 

Jy 


IC 797 


11 


58±0 


04 


11 


07±0 


06 


10 


80±0 


07 





200± 0.043 


0.730± 


045 


2 


110± 0.130 







085±0.021 


IC 800 


11 


19±0 


03 


11 


34±0 


03 


10 


58±0 


06 





130± 0.038 


0.600± 


056 


1 


560± 0.271 







076±0.019 


NGC 803 


10 


71±0 


02 


10 


01±0 


03 


9 


83±0 


04 





140± 0.040 


0.730± 


047 


2 


850± 0.134 


0.631±0. 22085 





093±0.019 


NGC 3270 


10 


76±0 


02 


10 


06±0 


01 


9 


73±0 


02 





130± 0.032 


0.620± 


040 


2 


280± 0.107 







059±0.014 


NGC 4712 


10 


92±0 


04 


10 


42±0 


05 


10 


09±0 


05 





070± 0.033 


0.500± 


044 


1 


880± 0.113 







102±0.023 


NGC 6137 


10 


51±0 


02 


9 


81±0 


02 


9 


55±0 


03 





020± 0.020 


O.lOOi 


100 





380± 0.380 







029±0.010 


NGC 7047 


11 


13±0 


03 


10 


47±0 


04 


10 


20±0 


05 





120± 0.025 


0.580± 


059 


2 


110± 0.251 







055±0.013 


NGC 7653 


10 


71±0 


02 


10 


03±0 


02 


9 


75±0 


04 





280± 0.041 


1.440± 


065 


5 


460± 0.434 







112±0.020 


PGC 035952 


12 


26±0 


05 


11 


58±0 


07 


11 


29±0 


07 





030± 0.034 


0.420± 


047 


1 


600± 0.212 


0.421±0. 14735 





051±0.013 


UGC 10205 


10 


87±0 


02 


10 


14±0 


02 


9 


89±0 


03 





060± 0.041 


0.460± 


056 


1 


560± 0.192 







058±0.015 


UGC 12519 


11 


62±0 


03 


10 


89±0 


03 


10 


61±0 


06 





130± 0.034 


0.960± 


052 


3 


280± 0.321 


0.408±0.1428 





074±0.016 



Note. — J, H a nd K magnitudes come from the 2MASS archive at |http: //irsa. ipac . caltech. edu/| the IRAS measurements were obtained using the 
Scanpi tool at |http : // scanpiops . ipac . caltech. edu: 9Q00/applications/Scanpi/index.html| The 450 and 850 fim measurements come from Vlahakis 
ct al. (2005) 
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must be subtracted prior to fitting the ellipses. 

For the IRAC data, the initial parameters were ob- 
tained from the Bertin & Arnouts (1995) source extrac- 
tor (SExtractor) catalog derived for each image, while the 
segmentation image produced by SExtractor was used to 
mask out stars and background galaxies. To produce a 
valid "data quality file" input to ellipse all pixels belong- 
ing to the object and the sky were set to zero in this im- 
age. In addition to masking foreground and background 
sources, the "data quality file" is also used to flag pixels 
affected by instrumental defects and cosmic rays. The 
same procedure to mask background/foreground objects 
was used for the MIPS observations. In the case of MIPS 
160 /im images, gaps in the sky that were not covered by 
the array were masked out by hand. 

Prior to the photometry, fluxes were converted from 
MJy/Sr into Janskys using the factors presented in Ta- 
ble 2, which depend on the pixel size. The sky level in 
the IRAC images was determined by fltting a low or- 
der polynomial in regions more than 100" away from the 
galaxies, which was then subtracted from the image. For 
MIPS images the background sky was estimated using 
several boxes of 25 pixels from which a mean sky value 
was determined and subtracted. For most galaxies the 
surface brightness proflle flattens at large radii, demon- 
strating that the sky subtraction was effective. However, 
in a few cases, in particular for 160 /xm, the presence 
of residual sky (caused by the foreground Galactic Cir- 
rus) can affect the final measurements of galaxies with 
very low fiux in this band. The circular and elliptical 
isophotes were fit to distances of ~ 100 " from the center 
and for most bands convergence was reached at about 
50". For most of the sources considered in this work, 
the results derived from elliptical and circular apertures 
differ by less than the estimated uncertainties. Based 
on this, the circular aperture photometry was adopted 
for the final flux measurements, which are presented in 
Table 3. 

2.3. Archival data of other wavelengths 

The Spitzer IRAC and MIPS data are supplemented in 
this paper with published data, presented in Table 4, con- 
sisting of the 2MASS photometry in the near- infrared, 
the sub-mm measurements by Vlahakis et al. (2005) and 
IRAS archival data. 

The 2MASS (Jarrett et al. 2000) near-Infrared mag- 
nitudes were obtained from the NASA Extragalactic 
Database (NED). The 2MASS measurements are the 
total magnitudes measured by extrapolating the near- 
IR surface brightness to a diameter corresponding to 
roughly 4 disk scale lengths defined from the J-band im- 
age due to its greater sensitivity at lower surface bright- 
nesses (Cutri et al. 2006). Because most of the galaxies 
presented in this work are relatively bright, the typical 
uncertainties in the photometry are ~ 3%. Because of 
the bright apparent magnitudes of these galaxies, the loss 
of light coming from low surface brightness regions must 
be comparable to that found by Bell et al. (2003) for 
K^13.5 galaxies, which is ~ 0.1 magnitudes, and should 
not affect our conclusions. 

The sub-mm data were presented by Vlahakis et al. 
(2005). The 850 /im fluxes were measured from the 
SCUBA maps choosing an aperture that would optimize 
the signal to noise ratio, taking into account the opti- 
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Fig. 2. — Spectral energy distribution for the ensemble of SLUGS 
galaxies observed for this project, normalized at 3.6 ^m. The stel- 
lar component shows little variation, while all galaxies, with the 
exception of NGC 6137 show a bump in the mid-IR due to the 
aromatic features, followed by a strong peak in the far-IR caused 
by the presence of cold dust. For NGC 6137 the Spitzer measure- 
ments at 70 and 160 ^m are only upper limits. The combined 24 
and 850 ^m data points are suggestive of a power-law emission in 
the mid- and far-IR for this galaxy. 

cal extent of the galaxies indicated by the Digitized Sky 
Survey*. 

The IRAS fluxes were obtained by co-adding images 
using the NASA/IPAC Infrared Science Archive's IRAS 
Scan Processing and Integration (SCANPI) tool and se- 
lecting the Noise- Weighted Mean. Given the much higher 
quality of the 24 /im fluxes from MIPS, we will not use 
the IRAS 25 /zm fluxes in the analyses that follow. 

3. ANALYSIS 

Figure 1 presents a mosaic of all images taken with 
IRAC and MIPS for the SLUGS galaxies observed with 
Spitzer. The images are ordered in observed wavelength 
and the scale bar in each panel represents a projected an- 
gular distance of 2 kpc. All galaxies are resolved in the 
IRAC bands, and different subcomponents of galaxies are 
easily identified. Similar to the SINGS galaxies (Dale et 
al. 2005) and the galaxies of Hinz et al. (2006, 2007), the 
older stellar populations provide the main contribution 
in bluer IRAC bands (3.6 and 4.5 /im) while the dustier 
disk becomes more prominent in the redder 5.8 and 8.0 
/tm bands. Although most galaxies of this sample are 
resolved in the MIPS bands, galaxy subcomponents are 
only easily identified in the MIPS 24 /im images. While 
the galaxy outlines and most prominent components can 
be seen in the 70/im images, the 160 /im observations in 
general show very little detail, though the extent to which 
the emission is seen is comparable to that of the stellar 
component. One galaxy (NGC 6137) is only weakly de- 
tected at 70 and 160 /im and shows a point-like source 
at 24 /im, probably due to the radio source at the center 
of this galaxy (Giovanini et al. 2005). The new data 
suggest that probably all, rather than ~ 20%, of the 850 

^ The Digitized Sky Surveys were produced at the Space Tele- 
scope Science Institute under U.S. Government grant NAGW-2166. 
The images of these surveys are based on photographic data ob- 
tained using the Oschin Schmidt Telescope on Palomar Mountain 
and the UK Schmidt Telescope. The plates were processed into 
the present compressed digital form with the permission of these 
institutions. 



7 



/xm emission detected by Vlaliakis et al. (2005) for this 
galaxy is due to an active nucleus. 

The spectral energy distribution of the SLUGS galax- 
ies normalized at 3.6 fxm, where the galaxy light is dom- 
inated by old stellar populations, is shown in Figure 2. 
Each galaxy is represented by a different color as noted in 
the key. Most objects show similar spectral shapes, with 
intensities peaking at ~ 100 /zm. The only exception is 
NGC 6137. 

3.1. Dust component temperatures 

The new MIPS data combined with existing IRAS and 
SCUBA measurements can be used to estimate the dust 
component temperatures and masses using composite 
grey-body curves (e.g., Vlahakis et al. 2005). For two 
components these are: 

F, = A^,y^B,{j^,T^) + A,iyl'B,{j^,T,) (1) 

where A,;, and Ac are the relative contributions due to 
the warm and cold dust components. The masses can be 
estimated from (Vlahakis et al. 2005): 



M = 



kd 



A,, 



Ar 



1 + 2' 



1 



Aw + An 



(2) 

where kd = 0.07 ±0.02 m^ kg ^ is the average dust mass 
opacity measured at 850 /xm by James et al. (2002), D 
is the galaxy distance and B^, the modified Planck func- 
tion. Since the average radiation field combines emission 
at several intensity levels (e.g., Draine & Li 2007; Draine 
et al. 2007; see below), the temperatures should only 
be taken as a rough estimate of the "effective dust tem- 
perature" (Blain et al. 2004). The fits were calculated 
considering three values for the dust cmissivity index (/3 
= 1.0, 1.5, 2.0), P = 2.0 producing better constrained 
fits. This is also the value of (3 used by Vlahakis et al. 
(2005) in their analysis. The fits are presented in Fig- 
ure 3, while Table 5 lists the temperatures and errors for 
the cold and warm components, and the estimated mass 
ratio between warm and cold components. The larger 
wavelength range of the present sample allows more con- 
strained fits than was possible for Vlahakis et al. (2005). 
For the 3 objects in common for which they calculate two- 
component fits, the most dramatic difference is found for 
UGC 12519, where the measurements change from (T^, 
Te, Mcoid/Mwarm) = ( 28, 18, 12) to (56.0, 19.6, 1054). 
For NGC 803 Vlahakis et al. (2005) find two possible 
solutions (33, 18, 92) and (60, 19, 2597). The latter 
is closer to the values measured here (54.7, 28.4, 1812). 
For PGC 35952 the IRAS+SCVBA measurements pro- 
duce (58, 18, 1859), while we find (54.7, 19.0,1054). The 
temperatures agree within the estimated errors, though 
the dust mass ratios can differ significantly. As noted by 
Draine et al. (2007) even small shifts in temperatures im- 
ply very large changes in the dust masses, as seen in this 
comparison. Because of the paucity of data longwards 
of 160 /xm for the SLUGS sample analysed here (only 3 
galaxies have 450 fim data) we cannot test for the exis- 
tence or absence of dust colder than 18 K which has been 
found in low-metallicity dwarf galaxies (e.g. Galliano et 
al. 2005). 



More robust estimates of the (total) dust masses can be 
obtained using models with a more complex dust compo- 
sition (e.g.. Desert et al. 1990; Guiderdoni et al. 1998; 
Dale et al. 2001; Dale & Helou 2002; Marshall et al. 
2006; Siebenmorgen & Kriigel 2007; Draine & Li, 2007; 
Draine et al. 2007). These models, combined with the 
empirical SEDs constructed using existing observations 
(e.g.. Chary & Elbaz 2001; Lagache et al. 2003; La- 
gache et al. 2004; Rieke et al. 2009) have been used 
to explain the strong evolution as a function of redshift 
in number and density that is observed in the mid and 
far infrared (e.g., Puget et al. 1995; Le Floc'h et al. 
2005; Perez- Gonzalez et al. 2008). Typically the mod- 
els describe the mid- and far-IR properties of galaxies 
combining the UV output from hot young stars, which is 
re-radiated by small size particles, and larger grains that 
are heated stochastically or are in thermal equilibrium 
(Desert, Boulanger & Puget 1990). 

To analyse the sample presented in this paper and 
compare their properties with those measured for SINGS 
galaxies, we follow the analysis of Draine et al. (2007) , by 
using the models calculated by Draine & Li (2007). The 
latter work builds upon many previous papers by Draine 
and collaborators (e.g., Li & Draine 2001; Weingartner & 
Draine 2001a; 2001b) by taking into account the results 
of recent observations (e.g., Smith et al. 2007). These 
models are based on a series of physically motivated as- 
sumptions, and can be used to estimate the galaxies' to- 
tal dust masses. The dust models are parameterized both 
by the dust composition and by the stellar radiation field 
incident on the dust grains. In the Draine et al. (2007) 
models, most of the dust is assumed to be located in a 
diffuse ISM that is heated by a general interstellar radia- 
tion field (e.g.. Lonsdale & Helou 1987) typical of normal 
spiral galaxies {Umin)- This field has an equilibrium tem- 
perature ~ 16.2 K (based on Milky Way measurements 
by Finkebeiner, Davis & Schlegel 1999) while a range of 
radiation field intensities heats the dust particles to dif- 
ferent temperatures. The fraction of galaxy dust that 
is subject to the intense radiation field generated by hot 
young stars (Jlmax) is denoted by 7 and is quantified as a 
percentage of the total dust emission of the galaxy. The 
emission at A < 20 /xm (Draine & Li 2007) is mainly due 
to aromatic grains, whose contribution is characterized 
by the qpah parameter, which quantifies the ratio be- 
tween the mass in aromatic grains containing less than 
1000 carbon atoms and the total dust mass. Weingart- 
ner & Draine (2001) calculated several models that use 
different grain size distributions and with different qpAH 
values, but are constrained to follow the average ISM ex- 
tinction law, reproducing the Milky Way's wavelength- 
dependent extinction curves (Draine & Li 2007). The 
models, described in Table 3 of Draine et al. (2007), 
show slight variations in the total dust mass to hydro- 
gen ratio but have values of Mdust/Mn ~ 0.01, slightly 
higher than the prefered value of Draine et al. (2007) 
of 0.0073, calculated assuming an approximately solar 
metallicity for the Galaxy's ISM. The stellar continuum 
for wavelengths greater than 5 /xm is well described by 
a 5000 K black-body (e.g.. Smith et al. 2007; Draine et 
al. 2007). Following Draine et al. (2007) the SEDs of 
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TABLE 5 

Temperature fit parameters for SLUGS and SINGS galaxies 



Galaxy 










K 


K 





IC 797 


53 


.9± 


1. 


,3 


18, 


.1± 


1 


.0 


1213 


IC 800 


51 


.1± 


2. 


,8 


16, 


,2± 


2, 


.8 


1065 


NGC 3270 


56 


.9± 


1. 


,9 


18, 


,5± 





,9 


1918 


NGC 4712 


56 


.2± 


0. 


,9 


17, 


.0± 


0, 


.4 


3213 


NGC 7047 


54 


.5± 


2. 


,8 


18, 


.8± 


1 


.4 


1248 


NGC 7653 


54 


.7± 


4. 


,0 


19, 


.9± 


2, 


.1 


1008 


NGC 803 


54 


.7± 


1. 


,0 


18, 


.4± 


0, 


.4 


1812 


PGC 035952 


54 


.7± 


4. 


,4 


19, 


.0± 


2, 


.5 


1054 


UGC 10205 


52 


.5± 


2. 


,8 


18, 


,4± 


1 


.8 


1098 


UGC 12519 


56 


.0± 


2. 


1 


19, 


.6± 


1 


.0 


1054 


NGC 337 


52 


.8± 


3. 


,2 


20, 


.6± 


3, 


.8 


337 


NGC 1482 


56 


.1± 


3. 


6 


23 


,5± 


5, 


.2 


109 


NGC 2798 


58 


.0± 


4, 


3 


24, 


,4± 


5 


,4 


116 


NGC 2976 


54 


.9± 


3. 


7 


19 


,1± 


5, 


.1 


800 


NGC 3190 


54 


.2± 


3. 


,4 


19, 


,9± 


3, 


.9 


766 


Mrk 33 


60 


.7± 


5. 


1 


24, 


.3± 


7, 


.2 


111 


NGC 3521 


55 


.2± 


2. 


,9 


18, 


,5± 


3, 


.8 


985 


NGC 3627 


55 


.2± 


3. 


,0 


19, 


.2± 


4, 


.6 


571 


NGC 4536 


55 


.5± 


2. 


,9 


20, 


.6± 


4, 


.6 


225 


NGC 4569 


58 


.4± 


3. 


,4 


19, 


.7± 


2, 


.9 


973 


NGC 4631 


52 


.6± 


3. 


,4 


18, 


.9± 


4, 


.0 


539 


NGC 4826 


52 


.4± 


3. 


,8 


20, 


.3± 


5, 


.0 


362 


NGC 5195 


59. 


4±27 


.9 


25, 


.7± 


9, 


.5 


258 


NGC 5713 


54 


.9± 


2. 


9 


20, 


.2± 


4, 


.1 


307 


NGC 5866 


57. 


2±16 


.1 


24, 


.5± 


3, 


.1 


878 


NGC 7331 


54, 


.4± 


3. 


,2 


19, 


.4± 


2, 


.8 


660 



TABLE 6 

Results using Draine & Li dust models fitting U, 



Galaxy Mdust Ldus* 7 <U> qpAH Umax 

M0 Lq % % 



IC 797 


6.71 


8.85 


7.11 


1.01 


4.58 


10^ 


IC 800 


7.51 


9.44 


15.55 


0.63 


3.90 


10^ 


NGC 3270 


8.33 


10.57 


0.55 


1.28 


4.58 


10" 


NGC 4712 


8.14 


10.17 


1.20 


0.79 


3.90 


10^ 


NGC 7047 


8.00 


10.37 


1.94 


1.73 


4.58 


10* 


NGC 7653 


8.05 


10.45 


4.55 


1.88 


4.58 


10=* 


NGC 803 


7.43 


9.64 


0.00 


1.20 


4.58 


10** 


PGC 035952 


7.77 


9.85 


8.57 


0.88 


2.50 


10* 


UGC 10205 


8.18 


10.38 


7.33 


1.16 


3.90 


10^ 


UGC 12519 


7.95 


10.31 


1.10 


1.71 


3.90 


10" 



TABLE 7 

Results using Draine Sz Li dust models with fixed Umax='^0' 



Galaxy 






7 


<U> 


gpAH 




M@ 


Lq 


% 




% 


IC 797 


6.57 


8.82 


0.74 


1.31 


4.58 


IC 800 


7.05 


9.40 


0.73 


1.64 


4.58 


NGC 3270 


8.33 


10.57 


0.55 


1.28 


4.58 


NGC 4712 


8.14 


10.17 


0.82 


0.78 


3.90 


NGC 7047 


8.02 


10.36 


0.65 


1.62 


4.58 


NGC 7653 


7.98 


10.44 


0.49 


2.12 


4.58 


NGC 803 


7.43 


9.64 


0.00 


1.20 


4.58 


PGC 035952 


7.51 


9.81 


1.86 


1.48 


2.50 


UGC 10205 


7.91 


10.36 


0.18 


2.04 


4.58 


UGC 12519 


7.95 


10.31 


1.10 


1.71 


3.90 
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Fig. 3. — Temperature fits using modified blackbody curves witii tiie emissivity index (3=2, assuming two temperature components. The 
dotted curves represent the individual components, the solid curve the combined flux. The values of best flt temperatures and associated 
errors are also shown in the figure. 
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SLUGS galaxies are fit by minimizing 



X 



E 



err' 



hs 



■ errt. 



(3) 



odel 



where, using the Draine et al. (2007) notation, the ob- 
served flux is a combination of the stellar component with 
the emission due to the diffuse ISM and dust heated by 
hot stars: 

Frnodel {vi) = fl*Biy{T^, Ui) + . '^^1 X 



47rL>2 

[(1 - j)pUjni,Urmn) + 



mini Umax-! '^)] 



(4) 

and pi, represent the power radiated by the dust mix- 
ture per unit mass and frequency that is subject to radi- 
ation field intensity of Umin or to a range of inten- 
sities {pv)^ assumed to follow a power-law characterized 
by the exponent a (Draine et al. 2007). Following these 
authors, we adopt a=2. The models for the PAH abun- 
dance, the radiation field, the total Mdust/Mn ratio and 
the extinction are combined into j„n the specific power 
radiated per unit dust mass^. 

When calculating the fits, we first normalize the stellar 
component (Q^) using the 2MASS photometry in J and 
H. This is followed by the calculation of the dust model 
normalization; final fits are obtained by varying 7 and 
looping through all models sharing the same proportion 
of aromatic masses (qpAH) and general radiation field 
(Umin)- The results of this procedure are presented in 
Table 6, while Figure 4 shows the measured SEDs and 
the best fitting model. Also shown in Fig. 4 are the 
graphic estimates of the stellar flux, the flux due to dust 
in photo-dissociation regions and the flux emitted by dust 
in "the diffuse ISM heated by a single radiation field" 
(Draine et al. 2007). Most galaxies are reasonably well fit 
by the cold temperature diffuse component with a small 
contribution from hot dust in photo-dissociation regions 
(PDR). In the case of NGC 803 the fits are consistent 
with no presence of PDRs, though the IRAC data do 
not show any gross inconsistency with the presence of 
aromatic emission, which is generally produced in regions 
with strong UV emission. The analysis by Bendo et al. 
(2008) of the correlation observed between the aromatic 
emission at 8 /xm and 160 /im in a sample of nearby 
spiral galaxies as well as that of Haas, Klaas & Bianchi 
(2002) shows that such emission is expected. No fits 
are presented for NGC6137 because its FIR spectrum 
is probably dominated by non-thermal emission. The 
results in Table 6 also fit for Umax, which represents the 
upper cut off of the stellar intensity field (Draine et al. 
2007). Fixing the value of Umax = 10^ produces the fits 
presented in Table 7. When fitting for Umax we find that 
the residuals in the far-IR are slightly smaller. 

3.3. Comparison with SINGS 

3.3.1. Sample biases 

The combination of selection criteria (in addition to 
the small sample sizes of SLUGS and SINGS), imposes a 
major limitation on the possible analyses with these sam- 
ples. As described by Vlahakis et al. (2005), the parent 



° Tables with jm values are available at 
[ftp : //ftp ■ astro ■ pr inceton . edu/dralne/d ust /ireiii4/ site. 



the 



sample for SLUGS is defined from the CfAl Redshift 
Survey (Huchra et al. 1983), which in principle should 
be unbiased as far as the cold-dust properties of galaxies 
are concerned (cf., Section 1). The SLUGS galaxies se- 
lected for Spitzer observations are a subset of this larger 
sample. 

The optical SLUGS sample was defined considering the 
number of objects within a well defined region of the sky 
and with an angular size that would fit inside the SCUBA 
field of view. The parent sample contains a total of 97 
galaxies, observations having been taken for 81. Of the 
latter, 52 were detected in 850 /Ltm, while 29 (36%) only 
provided upper limits. Early type (E-I-SO) galaxies com- 
prise 14 (48 %) of the non-detections (corresponding to 
17% of the total observed sample), spirals 8 (28%) and 
the remainder are either multiple systems or compact 
sources. The sample observed with Spitzer contains 10 
spiral galaxies and 1 elliptical, which is close to the mor- 
phological mix of the SLUGS observed sample. Thus the 
small proportion of early-types in the Spitzer sample 
reflects in part the low detection rate of these galaxies 
in the SLUGS sample (Vlahakis et al. 2005). There is 
a risk that the SLUGS / Spitzer sample is not represen- 
tative of cold-dust properties of the local population of 
galaxies, being biased towards galaxies with larger cold 
dust masses. 

In a similar way, SINGS galaxies by design were not 
chosen to represent a volume-limited sample of galaxies 
but rather to be representative of the local morphology, 
luminosity and Lpm/ Lb parameter space (Kennicutt et 
al. 2003). Thus, when interpreting the comparisons that 
follow one must take into account that neither sample is 
a "fair" representation of the local mix of galaxy proper- 
ties, so that the results may not be valid for the general 
galaxy distribution. In particular, the observations of 
SINGS galaxies with SCUBA did not follow any type of 
preselection based on their dust properties. Of the 26 
SINGS galaxies with SCUBA data, 13 were observed af- 
ter the SINGS sample was defined, while data for the 
remainder were taken from the SCUBA archive by Dale 
et al. (2005). However, we should note that Draine et 
al. (2007) identify 9 objects which have lower confidence 
SCUBA observations which are galaxies for which the 
data reduction removes extended emission (4 galaxies); 
galaxies for which SCUBA measurements do not cover 
the total spatial extent where emission in 160 /xm is de- 
tected (4 galaxies) and one galaxy for which the 850 /xm 
emission is dominated by an AGN. 

3.3.2. Comparison between SINGS and SIUGS 

One motivation for the new Spitzer observations in 
this work was the significantly different distribution ob- 
served using the sample of SLUGS galaxies relative to 
SINGS when considering the IRAS 25 to 60 /xm and 
100 to 850 /im color-color distribution, which measures 
the relative intensity of warm to cold dust emission. 
Figure 5 shows the analogous distribution using MIPS 
data for both samples in the 24 to 70 /im and 160 to 
850 /im colors, where it is clearly seen that SLUGS 
(solid circles) and SINGS galaxies with SCUBA mea- 
surements (diamonds) occupy mainly different regions 
(we should note that the parent SLUGS sample covers 
the entire IRAS 25 to 60 /im color range). Even when 
including SINGS galaxies for which the 850 /im may 
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Fig. 4. — Spectral energy distribution of SLUGS galaxies including data from 2MASS, IRAS, SCUBA, IRAC and MIPS observations 
(open circles). The (estimated) error bars for most measurements are smaller than data points for most of the IRAC and MIPS data. Colors 
represent fits using the Draine Sz Li (2007) models; the dashed orange line represents the estimated stellar flux, the green line the emission 
due to grains heated in photo-dissociation regions, cyan the emission from dust due to a single radiation field. The black line represents the 
sum of these components. The parameters for these fits, which also fit for Umax the maximum starlight intensity, are presented in Table 6. 
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Fig. 5. — Color-color distribution for galaxies with SCUBA 
measurements. SLUGS galaxies are represented by black solid 
circles while the diamonds denote SINGS galaxies with good 
SCUBA measurements. The arrows represent upper limits of the 
uFu(160fini)/uFu{850^im) ratio of galaxies for which the SCUBA 
measurements are probably underestimated (Draine et al. 2007). 
The segregation between both samples is very clear. While both 
sample show a similar behavior in the 24/70 fim color which mea- 
sures the mid-IR slope, the SLUGS galaxies show flatter 160/850 
fim slopes due to the larger contribution due to cold dust. 
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Fig. 6. — Color-color distribution following Draine et al. (2007) 
of SLUGS galaxies (black circles) compared to SINGS galaxies with 
good SCUBA (diamonds) and lower confidence (suares) detections, 
and without (triangles). The location of SLUGS galaxies towards 
lower values of (i'i^i/)70;im/(i'^i^)i60fim is expected for galaxies 
with significant amounts of cold dust (Draine et al. 2007). 

be underestimated, thus providing upper limits to the 
z^Fiy(160/im)/j/F^(850/im) ratio (noted by the arrows), 
the locus of SLUGS galaxies still differs from that of the 
increased sample. 

We now consider several combinations of measureable 
parameters to examine whether this difference found be- 
tween SLUGS and the SINGS subsample with SCUBA 
data also extends to full SINGS sample when galax- 
ies with no SCUBA data are included. The color-color 
plot using MIPS observations is presented in Figure 6, 
again showing SLUGS and SINGS galaxies. As in Fig. 
5, we discriminate galaxies with good SCUBA observa- 
tions (diamonds) from those with lower-confidence level 
measurements (squares) and no SCUBA observations 



(triangles). The plot shows that SLUGS galaxies are 
found towards lower values of i^Fi,(70/Ltm)/i/F;^(160/L(m) 
which is the expected behavior for galaxies with cold 
dust, though not reaching the levels of NGC 4725 and 
NGC 2841, which are the most extreme cases in the 
SINGS sample. On the other hand, the distribution of 
the i/Fi,(24/Ltm)/i/F^(70/zm) colors of the SLUGS sample 
is similar to that of the SINGS sample. 

The data in Table 7 are used to compare the proper- 
ties of SLUGS galaxies with the Draine et al. (2007) fits 
for SINGS where the maximum intensity of the starlight 
(Umax) is fixed to 10^. The SED fits calculated in Section 
3.2 suggest that most galaxies in the SLUGS sample have 
PAH emission, which is further corroborated by the anal- 
ysis of Bendo et al. (2008). The frequency distribution 
of the amount of dust contained in small aromatic grains 
{qpAn) is shown in Figure 7(a). The SLUGS sample 
is represented by the solid grey histogram while Draine 
et al. (2007) measurements for SINGS galaxies are out- 
lined as the dotted histogram for the entire sample and 
the solid line histogram for the SCUBA sub-sample of 
17 galaxies. The distribution of aromatic grain indices 
for SLUGS and the full SINGS is similar in spite of the 
small number of SLUGS objects. Figure 7(b) shows the 
distribution of 7, which quantifies the contribution due 
to starlight re-radiatcd in photo-dissociation regions rel- 
ative to the idealized single and constant radiation field 
Umin- The 7 distribution for SLUGS has a median value 
of ~ 0.6%, and the galaxies are all concentrated in a sin- 
gle bin in the histogram; SINGS galaxies have 7 ^ 1%. 
This low value of 7 for SLUGS has direct implications for 
the estimate of the average starlight intensity shown in 
Figure 7(c) which is lower (< U >Biugs^ 1-3) than what 
is measured for SINGS (< U >sings^ 2.5). 

Figure 8(a) shows the distribution of total dust masses. 
SLUGS galaxies have masses ranging from 10^'^ Af© 
to 10^-^ Mq, with a median of lO"^-^ Mq, ~ 3 x higher 
than is measured for the entire sample of SINGS galax- 
ies (10^'* Mq). However, as noted by Draine et al. 
(2007), for galaxies without sub- mm data the fitting pro- 
cedure has to be more constrained otherwise the result- 
ing fits may produce larger cold dust masses. How- 
ever, Draine et al. (2007) find that the mass estimates 
for ~ 65% of their sample agree within a factor of 1.5 
when fits using and discarding the sub-mm data are 
compared. Given the uncertainties, the observed differ- 
ence between SLUGS and SINGS is probably not signifi- 
cant, corroborated by the similar masses of SINGS galax- 
ies with sub-mm data (10^'^ Mq). Figure 8(b) shows 
the total dust luminosities, which for SLUGS galaxies 
{logioL/LQ ^ 10.2) are somewhat higher than for the 
entire SINGS sample [logi^L / Lq ^ 9.7) (though the 
distribution is virtually identical to that of the SINGS 
SCUBA subsample logioL/LQ ^ 10.3). Differences be- 
tween the distribution of the dust model parameters were 
examined using the two-sided Kolmolgorov-Smirnov test 
(e.g.. Press et al. 2002), and the only significant dif- 
ference between both samples is found for the average 
starlight field < U > where the probability of both dis- 
tributions being drawn from the same parent sample is 

On average, the SLUGS galaxies have a warm 
dust component with a median temperature of 
T,„=54.7±1.7 K and a colder one with Tr= 18.5±1.1 K. 
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qp,h(%) y (%) <u> 

Fig. 7.- Histograms comparing the parameters obtained by fitting the Draine et al. (2007) dust models to tlic SLUGS galaxies (solid 
grey histograms), with results tabulated by Draine et al. (2007) for the full SINGS sample (dotted line histogram) and the subsample of 17 
galaxies with SCUBA measurements (solid line histogram). This comparison uses the fit parameters obtained by fixing the value of Umax- 
Panel (a) shows the Qpah distribution, which measures the ratio between the dust mass in the form of aromatic grains and the total dust 
mass. Panel (b) shows the distribution of 7 which quantifies the contribution towards the total dust emission from dust located in regions 
with intense starlight due to very hot stars. Panel (c) shows the average starlight radiation field (<U>) present in the galaxies' interstellar 
medium. 
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Fig. 8. — Histograms comparing the estimated dust masses (left panel) and luminosities (right panel) for SLUGS (grey histograms), and 
SINGS galaxies (dotted line histogram). The 17 SINGS galaxies that have SCUBA measurements represented by the solid line histogram. 
The present sample has marginally higher dust masses than the SINGS sample, though the total dust luminosity is comparable. 



Both estimates agree within 2 <j with the results of Vla- 
hakis ct al. (2005) who find Ty,=A7A±2A K and Tc= 
20.2±0.5 K respectively. The warm dust temperature is 
comparable to that measured by Hinz et al. (2007) for 
UGC 6879, while the the average value we find for the 
cold dust component is comparable to the Hinz et al. 
(2006) measurement for the starbursting dwarf galaxy 
UGC 10445. Similarly to the findings of Vlahakis et al. 
(2005) and Hinz et al. (2007), most of the dust in these 
galaxies is cold, with Mcoidl M^arm ~ 1156±658. The 
distribution of the ratio between the cold and warm dust 
components for SLUGS and SINGS is shown in Figure 
9. Although there is some overlap in the mass ratios, 
SINGS galaxies with SCUBA measurements tend to have 
a smaller contribution due to cold dust with a ratio of 
538±305. 

The ratio between the total emission integrated from 
5-1000 /im, which includes the important contribution 
due to aromatic emission (e.g., Rieke et al. 2009) and the 
infrared luminosity in the range 100-1000 /um, which is 
mostly dominated by the cold dust component, is shown 
in Fig. 10. SLUGS galaxies are represented by the 



solid black circles, SINGS galaxies with SCUBA mea- 
surements by diamonds or open squares (in the case of 
low significance level measurements); the triangles rep- 
resent SINGS galaxies for which the far-IR is not as 
well constrained because of the absence of SCUBA data. 
Also shown are the sample of local LIRGs and ULIRGs 
(squares) and templates for lO''-^''^ < L < 10^^ Lq of 
Rieke et al. (2009) represented by stars. The SLUGS 
galax;ies populate a region of lower luminosities and with 
a higher proportion of far-IR emission than the sample 
of SINGS galaxies with SCUBA measurements. How- 
ever, it is interesting to note that the constrained Draine 
& Li (2007) models place many of the SINGS galaxies 
without SCUBA measurements in roughly the same pa- 
rameter space as the SLUGS galaxies. The combined 
SLUGS and SINGS and LIRG/ULIRG samples show a 
trend where the amount of energy in the far-IR becomes 
more important as lower luminosities are reached. This 
trend is also seen in the Rieke et al. (2009) templates. 

The new observations reported in this paper suggest 
that there is a wider variety of sub-mm properties than 
encompassed by SINGS alone, with a larger fraction 
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Fig. 9. — Distribution of the mass ratio between cold (T~19K) 
and warm (T ~ 55 K) dust estimated from modified blackbody fits 
for SINGS and SLUGS galaxies with SCUBA data. The filled his- 
togram represents SLUGS galaxies and the open histogram SINGS 
galaxies. Although there is significant overlap, SINGS galaxies 
tend to have smaller mass ratios compared to SLUGS. 
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(from the SLUGS selection). We combine new IRAC 
and MIPS data obtained with the Spitzer Space Tele- 
scope with existing 850 /xm, IRAS and near-IR mea- 
surements. The new measurements confirm the presence 
of cold (T'^ISK) dust in 10 of these galaxies. 

Using fits of Draine & Li (2007) models we show 
that the contribution due to light reprocessed by photo- 
dissociation regions is not as important in our sample as 
it is for SINGS galaxies. The models also suggest that 
the total IR emission of SLUGS galaxies has a small con- 
tribution due to PAH lines emitted in photo-dissociation 
regions. However the proportion of the total light con- 
tributed by the mid-IR emission is smaller than that 
measured for the SINGS sample. The typical total dust 
masses obtained using Draine & Li (2007) models are of 
lO^'^M© and dust luminosities are ~ 10 Lq. The most 
significant difference between the samples is the lower 
average intensity of the general radiation field of SLUGS 
compared to SINGS galaxies. Although the uncertainties 
are large, SLUGS galaxies show ~ 2 x larger proportion 
of cold dust relative to warm than do the SINGS galaxies. 

The present data combined with previous observations 
show that cold dust is present in all types of spirals be- 
ing the major contributor to their far-infrared luminosity 
(Tuffs and Popescu 2003). Although the 160 /im data are 
somewhat shallow, they suggest that qualitatively the 
scale length distribution of cold dust is at least compa- 
rable to what is measured in the mid-infrared and in the 
visible (e.g., Alton et al. 1998; Davies et al. 1999; Tuffs 
et al. 2002a, 2002b; Hinz et al. 2006). Since the amount 
of dust in a galaxy measures the quantity of heavy el- 
ements in the ISM, it gives a complementary measure 
of the average metallicity of galaxy (Dunne et al. 2000). 
Thus the presence of cold dust suggests that the extended 
HI in these galaxies is generally not metal-poor. 



Logio(L) (Lq) 

Fig. 10. — Ratio between the luminosity emitted in the far-IR 
(100-1000 ^lm) and the total luminosity (5-1000 fim) for SLUGS 
galaxies (black circles), SINGS galaxies with SCUBA detections 
(diamonds), those without (triangles); the Rieke et al. (2008) sam- 
ple of local LIRGs and ULIRGs (solid squares) and the Rieke et 
al. (2008) templates (stars) covering a range of luminosities from 
lO^ '^ to 10^3 Lq. The local LIRGs and ULIRGs extend to lower 
values the cold/total emission measured by SINGS and SLUGS. 
The behavior of the Rieke et al. (2008) templates show that the 
three samples are complementary. These also suggest that galax- 
ies with lower IR luminosities tend to have a larger contribution of 

f}Fga<ia«8§"itf t!lte feSifby universe presenting significant 
amounts of dust at temperatures ~ 18 K. The luminosi- 
ties of these galaxies (Log{L/LQ) ^ 10.2) are compara- 
ble to some of the least luminous objects detected in the 
SCUBA Half Degree Extragalactic Survey (SHADES), a 
blank field survey which scanned several patches of sky 
(Mortier et al. 2005), in contrast to SLUGS, which tar- 
geted optically-detected galaxies. The lower luminosity 
galaxies in SHADES are characterized as having "cirrus" 
SEDs (e.g., Clements et al. 2008) and make up 25 % 
of the sample of sub-mm galaxies (Clements et al. 2005), 
suggesting that the local samples are undersampling this 
type of galaxy. 

4. CONCLUSION 

We trace the Spectral Energy Distribution for 11 galax- 
ies providing an unbiased sample of sub-mm properties 
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